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The liquid-gas phase transition in hot neutron-rich nuclear matter is investigated 
within a self-consistent thermal model using different interactions with or without isospin 
and/or momentum dependence. The boundary of the phase- coexistence region is shown 
to be sensitive to the density dependence of the nuclear symmetry energy as well as the 
isospin and momentum dependence of the nuclear interaction. 



1. Introduction 

The liquid-gas (LG) phase transition in nuclear matter remains illusive and a 
hot research topic despite of the great efforts devoted to understanding its na- 
ture and exp erimental manifestations by the nuc lear p hysics community over many 
year a I I I I For a recent review, see, e.g., Refs! ^ ^ i Most of the previous stud- 
ies have focused on the LG phase transition in symmetric nuclear matter. While 
in an asymmetric nuclear matter, the LG phase transition is expected to display 
some distinctly new features because of the isospin degree of freedom and the as- 
sociated interactions and additional conservation laws . This expectation together 
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with the need to understand better properties of asymmetric nuclear matter rele- 
vant fo r both nuclear physics and astrophy sics have stimulated a lot of new work 
recentl^i^l^ ^^ l ^^ l ^-^ l ^^ l ^^ l ^^ l ^^ l ^^ l ^^ l ^ -^l^. Moreover, the study on the LG phase 
transition in asymmetric nuclear matter has received recently a strong boost from 
the impressive progress in developing more advanced radioactive beams that can 
be used to create transiently in terrestrial laboratories large volumes of highly 
asymmetric nuclear matter. Though significant progress has been made recently in 
studying properties of isospin asymmetric nuclear matter and the LG phase tran- 
sition in it, there are still many challenging questions to be answered. Among the 
main difficulties are our poor understanding about the isovect or nucle ar interaction 
and the density dependence of the nuclear symmetry energy I ^ ^ [ Fortunately, 
recent analyses of the isospin diffusion data in heavy-ion reactions have allowed 
us to put a stringent c onstraint on the symmetry energy of neutron-rich matter 
at sub-normal densities I ^4 | zo | /!d | j^ jg therefore interesting to investigate how the 
constrained symmetry energy may allow us to better understand the LG phase tran- 
sition in asymmetric nuclear matter. Moreover, both the isovector (i.e., the nuclear 
symmetry potential) and isoscalar parts of the single nucleon potential should be 
momentum dependent. However, effects of the momentum-dependent interactions 
on the LG phase transition in asymmetric nuclear matter were not thoroughly 
investigated previously. 

We report here our recent progress in investigating effects of the isospin and 
momentum dependent interactions on the LG phase transition in hot neutron- 
rich nuclear matter within a self-consistent thermal model using three different 
interactions^. The first one is the isospin and momentum dependent MDI interac- 
tion constrained by the isospin diffusion data in heavy-ion collisions. The second one 
is a momentum-independent interaction (MID) which leads to a fully momentum 
independent single nucleon potential, and the third one is an isoscalar momentum- 
dependent interaction (eMDYI) in which the isoscalar part of the single nucleon 
potential is momentum dependent but the isovector part of the single nucleon po- 
tential is momentum independent. We note that the MDI interaction is realistic, 
while the other two are only used as references in order to explore effects of the 
isospin and momentum dependence of the nuclear interaction. 



2. Theoretical models 

2.1. MDI interaction 

In the isospin and momentum-dependent MDI interaction, the potential energy 
density Vmdi(p, T, 5) of a thermally equilibrated asymmetric nuclear matter at total 
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density p, temperature T and isospin asymmetry 5 is expressed as follows l^^'^, 

Po 2po ^ cr + 1 pI 

PO^ 7 7 "^ l + (p-p')2/A2 

In the mean field approximation, Eq. JT]) leads to the following single particle po- 
tential for a nucleon with r noment um p and isospin t in the thermally equilibrated 
asymmetric nuclear matter I ° l ^ l 



Po Pa Pa 

X (1 - x5^) - 8t.t tL^Sp 



Po 7 " i + (p-f')2/a^ 

Po J l + {p-p'y/A^ ^ ^ 

In the above the isospin r is 1/2 for neutrons and —1/2 for protons, and fri^jP) is 
the phase space distribution function at coordinate r and momentum p. The detailed 
values of the parameters a, Au{x), Ai{x)^ B, Cr.r, Cr.-r and A can be found in Ref. 
| /a | zo | g^jj^^ have been assumed to be temperature independent here. The isospin and 
momentum-dependent MDI interaction gives the binding energy per nucleon of —16 
MeV, incompressibility Ki^ — 211 MeV and the symmetry energy of 31.6 MeV for 
cold symmetric nuclear matter at saturation density po = 0.16 fm~'^ ■ . The differ- 
ent X values in the MDI interaction are introduced to vary the density dependence of 
the nuclear symmetry energy while keeping other properties of the nuclear equation 
of state fixed l^. We note that the MDI interaction has been extensively used in 
the transport model for studying isospin effects in intermediate-energy heavy-ion 
collisions induced by neutron-rich nuclei ^-^ ^^ ^"^ "^^ '^^ '^'^ ^^ ^"^ ^^. In particular, 
the isospin diffusion data from NSCL/MSU have constrained the v alue of x to be 
between and —1 for nuclear matter densities less than about 1.2/9c|l "-* ! " I We will 
thus in the present work consider the two values of x = and x = —1 with x = 
giving a softer symmetry energy while x = —1 giving a stiffer symmetry energy. 
The potential part of the symmetry energy -Ef °^ (P' ^) ^^ ^^ro temperature can be 
parameterized byl^ 



El°l{p, x) = F{x)^ -f [18.6 - F{x)] (ii)«(-), (3) 

Po Po 



-^ -f [18.6 - i^(x)] (-^)^(") 
Po Po 

where the values of F{x) and G{x) for different x can be found in Ref. 
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2.2. MID interaction 



In the momentum-independent MID interaction, the potential energy density 
Vmid (Pi S) of a thermally equilibrated asymmetric nuclear matter at total density 
p and isospin asymmetry S can be written as 

Vmib{p, <5) = ^^ + -^^ + pEP^Up, x)5\ (4) 

2 Po 1 + 7 Po^ 

The parameters a, (3 and 7 are determined by the incompressibility Kq of cold 
symmetric nuclear matter at saturation density po as in Ref. ^ and Kq is again set 
to be 211 MeV as in the MDI interaction. And EP°l^(p,x) is just same as Eq. ^. 
So the MID interaction reproduces very well the EOS of isospin-asymmetric nuclear 
matter with the MDI interaction at zero temperature for both a; = and x — —1. 
The single nucleon potential in the MID interaction can be directly obtained as 

C/mid(p, 5, T)^a^+ P{^r + U^-^ip, S, r), (5) 

Po Po 

with 



U^^P,S,t) 



AF{x)-^ + 4(18.6 - F{x)){-^f^^^'^ 
Po Po 

(18.6 - F{x)){G{x) - 1)(-^)GW^2^ (6) 

PO 



2.3. eMDYI interaction 

The momentum-dependent part in the MDI interaction is also isospin dependent 
while the MID interaction is fully momentum independent. In order to see the 
effect of the momentum dependence of the isovector part of the single nucleon 
potential (nuclear symmetry potential), we can construct an isoscalar momentum- 
dependent interaction, called extended MDYI (eMDYI) interaction since it has the 
same f unc tional form as the well-known MDYI interaction for symmetric nuclear 
matterl^. In the eMDYI interaction, the potential energy density V^mdyi(p, T, 5) of 
a thermally equilibrated asymmetric nuclear matter at total density p, temperature 
T and isospin asymmetry 5 is written as 

A p^ B p^^" 

V'oMDYi(p, T, ^) — 



2 po 1 + o- po'" 

7oJ J ^^ l + ip-fy/A^ ^^sy™(P, ^)S ■ (7) 

Here fo{f,p} is the phase space distribution function of symmetric nuclear matter 
at total density p and temperature T. Again EPyl^{p,x) has the same expression 
as Eq. ©. We set A = ^"+^' and C = '^'^•-^^'^^■^ , eaidB, a and A have the same 
values as in the MDI interaction, so that the eMDYI interaction also gives the same 
EOS of asymmetric nuclear matter as the MDI interaction at zero temperature for 
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- 1 . The single nucleon potential in the eMDYI interaction can 
t/eMDYi (p, T, (5, p, r) = C/° (p, T, p) + C/-^^^ (p, 5, r ) , (8) 



where 



2C 



Po Po Po J l + ip-jfy/A^ 



h{f,p) 



(9) 



and U^^"^ (p, (5, r) is the same as Eq. ([5]) which implies that the symmetry potential is 
identical for the eMDYI and MID interactions. Therefore, in the eMDYI interaction, 
the isoscalar part of the single nucleon potential is momentum dependent but the 
nuclear symmetry potential is not. 

2.4. Thermodynamic Quantities 

At zero temperature, fT{r,p) — ■p-0(p/(T) — p) and all the integrals in above 
expressions can be calculated analytically'^, while at a finite temperature T, the 
phase space distribution function becomes the Fermi distribution 

2 1 



frif.p) 



/l3 



exp(- 



- + U^-l-i^ 



(10) 



1 



where pr is the chemical potential of proton or neutron and can be determined 
from 



fr{r,p)d'p. 



(11) 



In the above, m^ is the proton or neutron mass and Ur is the proton or neutron 
single nucleon potential in different interactions. From a self-consistency iteration 
scheme l*^"!^" !^ the chemical potential fir and the distribution function fT{r,p) can 
be determined numerically. 

From the chemical potential p,- and the distribution function f^ (r, p) , the energy 
per nucleon E{p, T, 5) can be obtained as 

„2 



E{p,T,5) 



V{p, T,S) + Y,J d'p£^Jrir,P) 



Furthermore, we can obtain the entropy per nucleon Sr{p,T^5) as 



Sr{p,T,5) 



Stt 



Ph'^ Jq 
with the occupation probability 



p [tIt Inrir + (1 — rir) ln(l — nr)]dp 



exp(- 



-+U^-Hr 



(12) 



(13) 



(14) 



) + l 



Finally, the pressure P{p, T, S) can be calculated from the thermodynamic relation 



P{p,T,6) 



Tj2Sr{p.T,S)-E{p,T,d) 



^PrPr- 



(15) 
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3. LG Phase Transition 

3.1. Chemical Potential Isobar 

With the above theoretical models, we can now study the LG phase transition in 
hot asymmetric nuclear matter. The phase coexistence is governed by the Gibbs 
conditions and for the asymmetric nuclear matter two-phase coexistence equations 
are 

Mf(T,p^,,S^)=Mf(T,p'^,J'^), (i = nandp) (16) 

P^(T,/,<5^)=pG(r,pG,<5«), (17) 

where L and G stand for liquid phase and gas phase, respectively. The chemical 
stability condition is given by 



dS 



> and ^ < 0. (18) 

The Gibbs conditions (|16p and (|17p for phase equilibrium require equal pressures 
and chemical potentials for two phases with different concentrations and asymme- 
tries. For a fixed pressure, the two solutions thus form the edges of a rectangle in the 
proton and neutron chemical potential isobars as a function of isospi n asym metry 
5 and can be found by means of the geometrical construction method II I 




0.0 0.2 0.4 0.6 0.8 1.0 
5 




Fig. 1. (Color online) The chemical potential isobar as a function of the isospin asymmetry 
5 at T = 10 MeV for the MDI and MID interactions(left) and the eMDYI interaction(right) 
with a; = and x = ~1. The geometrical construction used to obtain the isospin asym- 
metries and chemical potentials in the two coexisting phases is also shown. Taken from 
Ref.l2Zl. 



We calculate the chemical potential isobars at T = 10 MeV, which is a typical 
temperature of LG phase transition. The solid curves shown in the left panel of 
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Fig. [1] are the proton and neutron chemical potential isobars as a function of the 
isospin asymmetry 5 at a fixed temperature T = 10 MeV and pressure P = 0.090 
MeV/fm'^ by using the MDI and MID interactions with a; = and a; = — 1. The 
resulting rectangles from the geometrical construction are also shown by dotted 
lines in the left panel of Fig. [TJ When the pressure increases and approaches the 
critical pressure Pc , an inflection point will appear for proton and neutron chemical 
potential isobars. Above the critical pressure, the chemical potential of neutrons 
(protons) increases (decreases) monotonically with S and the chemical instability 
disappears. In the left panel of Fig. [1] we also show the chemical potential isobar 
at the critical pressure by the dashed curves. At the critical pressure, the rectangle 
is degenerated to a line vertical to the S axis as shown by dash-dotted lines. The 
values of the critical pressure are 0.265, 0.230, 0.195 and 0.154 MeV/fm'^ for the 
MDI interaction with x = 0, MID interaction with x = 0, MDI interaction with 
X = — 1 and MID interaction with x — —1, respectively. 

Shown in the right panel of Fig.[l]is the chemical potential isobar as a function of 
the isospin asymmetry 6 atT = 10 MeV by using the eMDYI interaction with x = 
and X = — 1. Compared with the results from the MDI and MID interactions, the 
main difference is that the left (and right) extrema of /in and /Zp do not correspond 
to the same 6 but they do for the MDI and MID interactions as shown in the 
left panel. The chemical potential of neutrons increases more rapidly with pressure 
than that of protons in this temperature. At lower pressures, for example, P = 0.090 
MeV/fm^ as shown in Panel (a), the rectangle can be accurately constructed and 
thus the Gibbs conditions P^ and ^T7\ have two solutions. Due to the asynchronous 
variation of /i„ and /ip with pressure, we will get a limiting pressure Pum above which 
no rectangle can be constructed and the coexistence equations (ITC)) and p?]) have 
no solution. Panel (b) shows the case at the limiting pressure with Piim — 0.205 
and 0.175 MeV/fm^ for x = and x — —1, respectively. With increasing pressure, 
in Panel (c) /i„ passes through an inflection point while ^j,p still has a chemically 
unstable region, and in Panel (d) /ip passes through an inflection point while /i„ 
increases monotonically with S. 



3.2. Binodal surface 

For each interaction, the two different values of S correspond to two different phases 
with different densities and the lower density phase (with larger S value) defines 
a gas phase while the higher density phase (with smaller S value) defines a liquid 
phase. Collecting all such pairs of 6{T,P) and 5'{T,P) thus forms the binodal 
surface. 

In Fig. [2] (a), we show the section of the binodal surface at T = 10 MeV for the 
MDI and MID interactions with x = and x = — 1. On the left side of the binodal 
surface there only exists a liquid phase and on the right side only a gas phase exists. 
In the region of "filet mignon" is the coexistence phase of liquid phase and gas phase. 
Interestingly, we can see from Fig. [2] (a) that the stiffer symmetry energy (x = — 1) 
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Fig. 2. (Color online) (a) The section of binodal surface at T = 10 MeV in the MDI and 
MID interactions with x = and x = —1. The critical point (CP), the points of equal 
concentration (EC) and maximal asymmetry (MA) are also indicated, (b) The section of 

binodal surface at T = 10 MeV in the eMDYl interaction with a; = and x = —1. LP 

27 
represents the limiting pressure. Taken from Ref. ' — '-. 



significantly lowers the critical point (CP) and makes the maximal asymmetry 
(MA) a little smaller. Meanwhile, the momentum dependence in the interaction 
(MDI) lifts the CP in a larger amount, while it seems to have no effects on the 
MA point. In addition, just as expected, the value of x does not affect the equal 
concentration (EC) point while the momentum dependence lifts it slightly (by about 
0.005 MeV/fm^). These features clearly indicate that the critical pressure and the 
area of phase-coexistence region in hot asymmetric nuclear matter is very sensitive 
to the stiffness of the symmetry energy with a softer symmetry energy giving a 
higher critical pressure and a larger area of phase-coexistence region. Meanwhile, 
the critical pressure and the area of phase-coexistence region are also sensitive to the 
momentum dependence. The MDI interaction has a larger area of phase coexistence 
region and a larger value of the critical pressure, compared to the result of MID 
interaction in the temperature of T = 10 MeV. 

Fig. [2] (b) displays the section of the binodal surface at T = 10 MeV by using 
the eMDYl interaction with x = and x = —I. We can see that the curve is cut 
off at the limiting pressure with Pu^ ~ 0.205 and 0.175 MeV/fm'^ for a; = and 
a; = — 1, respectively. We can also see that the limiting pressure and the area of 
phase-coexistence region are still sensitive to the stiffness of the symmetry energy 
with a softer symmetry energy {x = 0) giving a higher limit pressure and a larger 
area of phase-coexistence region in this temperature. 

Comparing the results of the MDI and MID interactions shown in Fig. [2] (a). 
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we can see that for pressures lower than the hmiting pressure, the binodal surface 
from the eMDYI interaction is similar to that from the MDI interaction. This 
feature implies that the momentum dependence of the symmetry potential has 
little influence on the LG phase transition in hot asymmetric nuclear matter while 
the momentum dependence of the isoscalar single nucleon potential significantly 
changes the area of phase-coexistence region for pressures lower than the limiting 
pressure. For pressures above the limiting pressure, the momentum dependence of 
both the isoscalar and isovector single nucleon potentials becomes important. 

4. Summary 

In summary, we have studied the liquid-gas phase transition in hot neutron-rich 
nuclear matter within a self-consistent thermal model using three different nu- 
clear effective interactions, namely, the isospin and momentum dependent MDI 
interaction constrained by the isospin diffusion data in heavy-ion collisions, the 
momentum-independent MID interaction, and the isoscalar momentum-dependent 
eMDYI interaction. At zero temperature, the above three interactions give the 
same EOS for asymmetric nuclear matter. By analyzing liquid-gas phase transition 
in hot neutron-rich nuclear matter with the above three interactions, we find that 
the boundary of the phase-coexistence region is very sensitive to the density depen- 
dence of the nuclear symmetry energy. A softer symmetry energy leads to a higher 
critical pressure and a larger area of the phase-coexistence region. In addition, the 
area of phase-coexistence region are also seen to be sensitive to the isospin and 
momentum dependence of the nuclear interaction. For the isoscalar momentum- 
dependent eMDYI interaction, a limiting pressure above which the liquid-gas phase 
transition cannot take place has been found. 
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